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The substructure of jets produced in an exclusive and a charm-induced dijet 
sample in photoproduction and in charged and neutral current interactions 
has been studied with the ZEUS detector at HERA. Jets were identified using 
the longitudinally invariant kj- cluster algorithm in the laboratory frame. The 
substructure of jets has been studied in terms of the jet shape and subjet 
multiplicity. Comparisons between the dijet sample and the quark-induced 
samples allow an extraction of the jet substructure for gluons. Leading-logarithm 
parton-shower Monte Carlo calculations give a good description of the differ- 
ences between quark- and gluon-initiated jets. In neutral current interactions, 
the measurements have been compared to next-to-leading-order QCD calcula- 
tions which are used to make a determination of the strong coupling constant, a 3 . 

1 Introduction 

The internal structure of jets gives insight into the transition between partons 
produced in the hard scattering process and the experimentally observable 
jets of hadrons. The substructure of jets is expected to depend mainly on 
the type of primary parton, quark or gluon, and to a lesser extent on the 
particular hard scattering process. QCD predicts that at high transverse 
energy of the jet, E^ 1 , where fragmentation effects become negligible, the jet 
structure is driven by gluon emission from the primary parton and gluon jets 
are broader than quark jets due to the larger colour charge of the gluon. The 
substructure of jets has been studied in terms of the integrated jet shape u, 
^(r), and the differential jet shape EJ, p(r) —where the energy flow inside a jet 
is considered, and the subjet multiplicity 13, (n s bj), where jet-like structures 
(subjets) within a given jet are studied. The data samples used in these 
analyses were collected with the ZEUS detector at HERA. During 1995-1997 
(1999-2000) HERA operated with positrons of energy E e = 27.5 GeV colliding 
with protons of energy E p = 820 GeV (E p = 920 GeV). 

2 Jet substructure in exclusive and charm-induced dijets in jp 

Jet shapes and subjet multiplicities in exclusive and charm-induced dijet sam- 
ples in photoproduction (PHP) for the kinematic region Q 2 < 1 GeV 2 and 
inelasticity 0.2 < y < 0.85 have been measured. Dijets events were iden- 
tified with the /cT-cluster algorithm u and selected with E^ 1-2 > 7,6 GeV 
and — 1 < rf et < 2. Charm quarks were tagged by identifying D*^ mesons 
through the K2ir decay mode using the AM method . The jet with closest 
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distance in azimuthal angle to the D* meson was associated with the charm 
meson. The jet not associated with the D*, referred to as the "untagged- 
charm" jet, represents the unbiased, i.e. not influenced by the D* selection, 
jet candidate whose internal properties are studied. 

The measured (n s bj) as a function of the resolution scale, y cu t, for the 
exclusive dijet sample is shown in Fig. 0a. The predicted (n s bj) is larger for 
gluon-initiated jets than for quark-initiated jets and the measured data are 
located between the two curves, -showing that the dijet sample is a mixture of 
quark and gluon jets. PYTHIA □ gives a good description of the data. Fig. |l|b 
shows (n s bj) as a function of y cu t for the "untagged-charm" jet. The agreement 
between data and theory is very good and the predictions of charm-initiated 
jets are consistent with the measurements. 

The dependence of the substructure of jets with rf et has also been studied. 
In the exclusive dijet sample, (n s f,j) increases with rf et (Fig. ||a), which is 
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Figure 2. (a) {n B hj) a * a fixed y C ut = 0.01 as a function of r) 3et corrected to hadron level 
for the exclusive dijet sample and the "untagged-charm" jet . (b) Integrated jet shape of 
the measured "untagged-charm" jet and the extracted gluon jet substructure. 
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Figure 3. (a) Comparison of the subjet multiplicity as a function of y C ut corrected to 
hadron level for the "untagged-charm" jet PHP sample and inclusive jet production in NC 
DIS . (b) Comparison of (n 3b j) for a fixed y cu t = 0.01 as a function of Ei e for inclusive jet 
production in NC and CC DIS. (c) Dependence of (n sb j) with rp^ 1 for a fixed y C ut = 0.01 
in CC DIS. 

consistent with the predicted increase in the fraction of gluon-initiated jets 
with rj° et . In the "untagged-charm" jet sample, the results are consistent with 
a pure sample of quark jets for —1 < rf et < 1.5. For the highest if et values, 
the data show a deviation from the prediction for quark-induced jets. Since 
the estimated gluon contamination to the charm-induced jet sample (mainly 
due to "charm excitation") has its highest contribution in the forward region, 
the deviation could be explained by the increase of the gluon-jet fraction in 
the charm-enriched sample. 

Since the dijet PHP sample consists of a mixture of quark and gluon 
jets, any measured observable O of the internal structure can be written as: 
O dijet = fq-O quar k + fg-Ogi uon , where fq, fg = l~f q arc the fractions of quark 
and gluon jets. The measurements of the substructure of the charm-enriched 
sample at high transverse energies (Ei^ > 15 GeV) can be considered as 
measurements for a pure sample of quark jets {O c ) larm — O quar k)- Taking the 
fractions f q , f g from LO Monte Carlo, the substructure of gluon jets can be 
extracted. The fraction f q predicted for a dijet sample with-|£^ e * > 15 GeV 
and -1 < r] jet < 2 by PYTHIA (/, = 0.66) and HERWIG U(f q = 0.69) are 
found to be similar. Fig. ||b shows the extracted gluon jet substructure, which 
is consistent with the QCD predictions. 

3 Jet substructure in deep inelastic scattering 

Quark initiated jets are expected to be predominant in charged current (CC) 
and neutral current (NC) deep inelastic scattering (DIS). In Fig. ||a, (n s bj) as 
a function of y cu t in the "untagged-charm jet" PHP sample and an inclusive 
jet NC sample (Q 2 > 125 GeV 2 ) are compared. The measurements are found 
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Figure 4. Measured (a) (n s ),j) and (b) (*I/(r)) for a fixed y C ut 
The data are compared to NLO QCD calculations for different values of c?g. 

to be similar and in very good agreement for resolution scales y cu t > 0.01, 
where the charm mass effects are negligible. The charm-initiated jets are very 
similar to light quark jets. Fig. [sJd shows the evolution of {n s i,j ) with E J T et for 
a fixed y cut — 0.01 in NC and CC processes. The value of (n s bj) decreases as 
E^ 1 increases. The agreement between both measurements indicates that the 
pattern of parton radiation within quark jets is to a large extent independent 
of the hard scattering process. In Fig. ||c (n s bj) for a fixed y cu t = 0.01 as a 
function of V jet rf^ CC interactions (Q 2 p > 200 GeV 2 ) is shown. Both ARI- 
ADNE (CDM) □ and LEPTO (MEPS) Q give a good description of the data. 
The substructure in DIS processes shows no dependence with "tf et . 

In NC processes, the subjet multiplicity (Fig. |^a) as a function of E^ 1 
for a fixed y cu t — 0.01 and the mean integrated jet shape (Fig. ^b) for a fixed 
radius of r = 0.5 have been compared to NLO QCD predictions and as has 



been determined. The extracted value from 



\ n sbj , 



for jets with E J T > 25 



GeV is a s {M z ) = 0.1185 ± 0.0016(stat)l^[!g(e^)to:uo7i(^)- 
The extracted value from (^(r)) for jets with E 3 ^ > 21 GeV is 

a s {Mz) = 0.1179 ± 0.0017( S tat)l° ;^(ea ;p )irou 9 73(^)- 

Both values are in good agreement with the world average value u. 
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